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PREFACE 


Most  researchers  have  used  finite  differences  with  nonlinear  scaling  to  solve 
the  problem  of  high  current  Injection  Into  bipolar  devices.  This  report  records 
an  attempt  to  solve  the  same  problems  with  another  technique  and  resolve  the 
problems  with  fewer  restrictions  than  have  been  used  previously. 

Chris  Ashley,  Roe  J.  Maier  and  A.  Brent  White  of  the  Electronics  Division, 

Air  Force  Weapons  Laboratory  (AFWL),  have  provided  much  useful  information 
necessary  for  the  work  reported  in  this  report.  Dr.  Donald  C.  Wunsch  (AFWL/EL) 
provided  invaluable  assistance  in  the  research  efforts  of  this  problem.  I wish 
to  thank  these  individuals  for  their  help. 
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SECTION  I 
INTRODUCTION 

This  technical  report  documents  the  results  of  an  effort  to  model  the  effects 
of  high  current  injection  into  bipolar  semiconductors.  The  work  was  started  at 
the  recommendation  of  Dr.  Donald  C.  Wunsch  of  the  Air  Force  Weapons  Laboratory. 

When  a transistor  device  experiences  an  electromagnetic  pulse  (EMP),  high 
currents  and  voltages  can  be  induced  between  the  metal  contacts  of  the  device. 
Depending  on  the  materials  and  doping,  significant  changes  can  occur  in  the  dis- 
tributions of  the  electric  field,  the  electron  density,  the  "hole”  density,  and 
traps  (ref.  1).  Because  of  the  significant  heating  of  the  devices  at  high 
current,  a temperature  dependency  was  included  in  the  model.  The  number  of 
spatial  dimensions  was  reduced  to  one  to  keep  the  problem  tractable. 

Two  techniques  or  models  were  used  in  an  attempt  to  solve  the  above  problem. 
The  first  technique  employed  finite  differencing  in  both  variables  (time  and 
space)  and  the  second  was  based  on  reducing  the  model  equations  to  a system  of 
first  order  differential  equations  and  then  solving  this  system  of  equations. 
Difficulties  were  encountered  in  each  of  the  two  techniques.  In  both  cases,  the 
difficulties  were  related  to  the  physical  problem.  This  problem  is  characterized 
by  very  large  boundary  values  which  affect  the  solutions  obtained.  Because  of 
this  problem  and  others  described  in  this  report,  no  successful  predictions  were 
made  describing  the  effects  of  high  current  injection. 
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SECTION  II 

DERIVATION  OF  THE  MODEL  EQUATIONS 

Various  assumptions  are  made  concerning  the  special  dimensions,  the  form  of 
the  coefficients,  and  the  number  of  dependent  and  independent  variables.  The 
first  assumption  is  that  of  one  special  dimension,  say  x.  Using  this  assumption. 


one  reduces  Maxwell's  equations  as  follows  (ref.  2): 

curl  E * - ll  (la) 

curl  H - 1^  + J (lb) 

div  B ■ 0 (Ic) 

div  D ■ p (Id) 

curl  E - |i  (0,  -1,  +1)  » - II  (Ic) 


where 


D(x,  t) 


» e(t  - f)  E(X,  f)  df 


div  B 


3B 

dx 


- 0 


from  the  dimensional  assumption,  then  one  may  conclude  that’curVH’® ’0~Th^jneH^ 
assumption  is  about  the  form  of  the  current  density  J.  One  usually  assumes  that 
the  current  density  takes  the  form  (refs.  3,  4,  5,  and  6) 

Jn  “ n + q Dn  H » q Uj  U^  U‘  + q Dj  1^  (2a) 
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Jp"'Vp.qOp|E=qM.  U>U‘-qD.|l 


(2b) 


with 


J = Jn  " Jp 

The  quantities  pn,  yp,  n,  p,  Dn,  Dp,  q are  defined  in  the  glossary  of  terms; 
through  U®  are  a relabeling  of  these  variables.  Some  constants  will  be  relabeled 
in  view  of  their  connection  with  the  various  dependent  variables. 

Since  curl  H = 0,  equation  (lb)  reduces  to. _the_  equation  _ _ __ 


3D 

It 


q D3  - q(y3  U*  U«  + y^  U«)  - q D, 


'2  3x 


and  the  generation  rates  for  and  are  given  by  the  equations  (refs.  3,  4, 
5,  and  6) 


3U* 

3t 


3U’ 

3t 


n . 1 3 . 


n U*  1 3 n ^ 
^ ■ T3  ■ q 3x  ''^p' 


(3a) 

(3b) 


where  g * g(x,  t,  u),  u * u (Uj,  Uj,  U3,  u^,  U5,  Ug).  The  terms  are  related 
to  the  trapping  of  electrons  and  holes  by  the  media.  The  are  simply  the 
average  life  times  of  holes  and  electrons  between  capture.  Substituting  equa- 
tions (2a)  and  (2b)  into  equations  (3a)  and  (3b)  results  in  the  following: 
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I- 

f 
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1^  “ - < Vj  aj  > U*  U*  + < Vj  aj  > (N  - U®)  (4e) 

U‘  = eU«  (4f) 

The  equation  u'  = eU*  is  an  approximation  that  is  valid  at  low  frequencies  and 
isotropic  media.  The  size  of  a large  bipolar  transistor  is  about  one  millimeter 
in  its  largest  dimension.  For  such  a transistor  to  "see"  or  sense  wave  phenomena, 
the  frequency  of  the  wave  must  be  larger  than  10"  hertz  since 

f ■ C/r  » (3  X 10*  m/sec)/10''*  m) 

■ 3 X 10"  hertz. 

This  is  in  the  lower  range  of  the  region  of  visible  light  (near  the  upper  infra- 
red region).  This  approximation  allows  one  to  reduce  the  number  of  variables 
from  six  to  five.  The  reduced  set  of  equations  is  as  follows: 


_ 3U‘ 

3t  ■ 

-q(y2  + 

U3  u*) 

u‘ 

+ ql 

au* 

ax 

. D 

2 ax  / 

(5a) 

11 

g(x,  t,  u) 

-yi. 

^2 

a 

ax 

U* 

u‘  + 

D |yi\ 

2 ax  J 

{5b) 

II 

g(x,  t,  u) 

U* 

" T,  ■ 

a 

ax 

(“• 

U* 

u*  - 

D 

ax  / 

(5c) 

3Ul  3 
at 

K a^u"  _ 
a aP 

q(U‘)" 

(Uz 

u' 

f u 

, u^) 

- qU‘  (dj 

au® 

ax 

■ D 

ax  / 

(5d) 

au® 

at 

- < V2  Oj 

> u® 

- 4 

• *1 

> u® 

U®  + < V, 

o,  > 

N u® 

(5e) 

The  terms  < V.  > are  the  average  cross  sections  for  the  electrons  (j  = 2)  and 
J J 

holes  (j  = 3),  i.e.,  capture  possibilities.  The  0^  are  diffusion  parameters  which 
depend  on  the  electric  field  u',  temperature  u"*,  and  the  two  concentrations. 
Equation  (5e)  controls  the  number  of  traps  and  number  of  electrons  available  to 
the  system  as  a whole,  i.e.,  it  gives  the  system  charge  conservation. 
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SECTION  III 

FINITE  DIFFERENCE  TECHNIQUE 

Finite  differences  was  the  first  of  two  techniques  by  which  the  solution 
of  equations  (5)  was  effected.  This  technique  makes  use  of  the  definition  of 
the  partial  derivative  of  a function,  i.  e. 


and 


and 


3f  /„  4.\  _ Lim  f(x,  t h)  - f(x,  t) 
3t  ^ “ h^  h 

~ f(x.  ^ ~ 


3f  _ Lim  f(x  + h,  t)  - f(x.  t) 
3x  ” h^O  h 


3*f  _ Lim  3x 


(x  + h,  t)  - 4 (x.  t) 


3x' 


h->0 


3*f 


A similar  expression  holds  for  by  partitioning  the  x-t  space,  one  can  step 
through  the  time  portion  of  the  axes  or  the  spatial  axis.  The  second  derivative 
of  f with  respect  to  (say  x)  is  replaced  by 

3x* 


(Ax)^ 


where 


'j.k  ■ •'(’'j-  ‘ki 


X + jAx 


t|^  ■ t + kAt 


k,j  « 1, 


M 


This  leads  to  the  following  set  of  explicity  equations  (note  that  superscripts 
refer  to  the  name  of  the  variable  and  subscripts  tell  where  in  the  mesh  the 
function  is  evaluated). 
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■ “j.k  ’ Ir  ■ (“j.k  “j,k  * “j.k  '’j.k)'^ 


* “j.k  (“j.k  - “j-l.k)  - “j.k  (“j.k  - “j-l.k) 

“j.k*l  - “j.k  ■ “ k («;_,,  - 

* “j.k  f “j.k  (“j.k  -“j-l.k)*  “j.k  (“j.k  -“j-l.k) 

* “j.k  ^(“jkl.k-  ““j.k*  “j-l.k) 

* ^ (“j.k  - “j-l.k)  (“j.k  - “j-l.k) 


“j.k*l  - “j.k  * “ 9j,k  - - ^'ix  “j.k  - “j-l.k 


^ (“j.k  - “j-l.k)  “j.k  - “’L  (“j.k  - “j-l,k) 


* “j.k  i7(“jkl.k-^  “j.k*  “j-l.k) 

*Z7  (“j.k -“j-l.k)  (“J.k -“j-l.k) 


“j.k+l 


“j.k^l 


_ K At 

U • I,  ” — — 


T (“j*l.k  - 2 “j.k  ^ “j-l.k)  * “j.k  (“j.k  - “j.k-i)““ 


- “j.k  Uj.k  “j.k  “j.k  “j.k 

♦ < ,,  c.  > N 
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which  is  the  explicit  formulation  of  equations  (5)  and  where 

9,-  ^ = gfxi  U)  K,  . ,,U?  K,  U!|  U?  I,,  U! 

V Jj  J>K  J**'  J»K  J»K/ 

" 9(^j’  ^k*  ^j.k) 


i.e.,  Uji^  = ujk.  By  explicit,  one  means  that  the  right  hand  side  of  the  equation 
depends  on  known  values  at  previous  mesh  points.  Richtmeyer  and  Morton  (ref.  7) 
have  shown  that  under  certain  circumstances  the  explicit  formulation  of  finite 
differences  will  be  unstable.  It  is  possible  to  rewrite  equations  (6)  so  that 
the  system  is  either  implicit  or  a combination  of  implicit  and  explicit  as  shown 
for  the  diffusion  equation  in  reference  8 by  introducing  a parameter  0.  For 
example,  equation  (6d)  becomes 


- u!: , = 0 


"j.k+1  '"j,k 


KAt 

Ax^ 


(’^j+l,k+r  ^ '^j.k+l  ^ '^j-l,k+l, 


where  0e[O,  1]9R.  For  0=0  the  above  explicit  formulation  is  obtained  while  for 
0 = 1,  a fully  implicit  method  is  obtained.  The  Crank-Nicholson  method  is 
obtained  when  0 = 1/2.  Ford  and  others  have  shown  that  the  truncation  error  for 
the  diffusion  equation  using  constant  coefficients  and  the  Crank-Nicholson  method 
is  0 {(Ax)^  + (At)^}  while  the  purely  implicit  technique  and  explicit  techniques 
give  0 {(ax)^  + (At)}  (refs.  7 and  8).  In  reference  7,  the  authors  show  that  for 
constant  coefficients  the  finite  difference  scheme  is  stable  for  0e  (1/2,  1). 

Lees  (ref.  9)  has  shown  that  the  Crank-Nicholson  scheme  is  stable  for  certain 
types  of  quasi-1 inear  parabolic  equations.  Lees  states  that  how  one  differences 
the  lower  order  terms  in  the  equation 

= Ai(x,t)  11+  A^(x,  t)  u(x,  t) 

8x 

is  a matter  of  truncation  error  and  not  stability.  There  are  requirements  on  a 
and  u that  are  needed  so  that  a pure  implicit  finite  difference  operator  will  be 
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stable  ^ order  to  determine  whether  or  not  the 

equations  (6)  are  stable,  the  coefficients  of  the  equations  of  second  order  will 
have  to  be  analyzed  to  determine  if  they  satisfy  the  criteria  given  by  Lees. 

The  first  order  equations  will  be  inspected  at  a later  time. 


The  parameters  u^  and  are  the  mobilities  and  diffusion  rates  of  the  elec- 
trons and  holes  in  the  medium.  Various  forms  of  equations  are  used  according  to 
some  of  the  hypotheses  made  by  the  author  of  the  paper.  Maier  (ref.  1)  reports 
on  using  a multitrap  model  as  well  as  one  trap  model.  Van  Lint  (ref.  6); 

Raburn  and  Causey  (ref.  3);  Leadon  (ref.  4);  and  Newdeck  (ref.  5)  all  use  a 
variation  of  the  model  due  to  Shockley  and  Read  (ref.  10).  These  equations  take 
the  form 


= _J[ 


a.exp  (-g./UM 


L 


(U“) 


(8a) 


Oj  = Vj/U‘  (or  Dj  = constant) 


(8b) 


depending  on  the  author.  In  many  cases  the  constraints  or  parameters  are  curve 
fitted  between  experimental  data  assuming  a particular  form  of  vai^iation.  The 
function  g is  the  recombination-generation  term  which  gives  the  relationship 
between  electron-hole  pairs  during  the  generation  and  recombination  processes. 

Using  the  Shockley-Read  model  (ref.  10),  the  form  of  g is  (refs.  3,  4,  and  5) 


g(x,  t,  U)  = 


T2  (U"  + N^)  + Tj  (U'  + N^) 


(9) 


wi  th 


= 0(0**)^/^  exp  (-q(ff  - E')/{2  KijU**)) 
and  a constant. 


The  equations  of  interest  are  (6b),  (6c),  and  (6d)  or  (4b),  (4c),  and  (4d).  The 
coefficients  of  the  lower  order  terms  in  equation  (6b)  are  y,,  y®  and  g - — + 

u lui  D 

y2  i “2 • 

•^x 
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The  requirement  given  by  Lees  (ref.  9)  for  the  terms  Aj,  Aj  is  that 


but  note  that  these  A,,  do  not  depend  on  the  solutions;  hence  we  may  not  be  able 

to  use  this  theory  unless  we  assume  that  the  terms  , g and  U*  are  all 

uniformly  bounded  or  constant.  In  reference  9,  Lees  requires  that  the  coeffi- 
^211 

cient  of  the  term  have  uniformly  bounded  partial  derivative  as  well  as  being 

oX^ 

uniformly  bounded  away  from  zero.  If  these  conditions  are  satisfied  then  the 
equation 

du  / *.  \ 3^u 

It  = 


will  be  unconditionally  stable  with  respect  to  the  norm  |I*|1* 
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SECTION  IV 

ORDINARY  DIFFERENTIAL  EQUATION  APPROACH 

Another  technique  for  solving  systems  of  partial  differential  equations  is 
discussed  by  Shampine  (ref.  12)  for  the  diffusion  equation.  Basically,  one 
finite  differences  the  special  derivatives  and  steps  through  in  time  or  finite 
differences  in  the  time  and  solves  Poisson's  equation.  The  choice  depends  on 
how  complicated  the  derivatives  are  and  the  form  of  the  resulting  matrices. 
Assuming  that  one  finite  differences  the  spatial  variables  then  equations  (4) 
become 

dul 

• q(i..ujuj  * u,Uj‘  uj)  * ,(oj(Uj>-  uj.,)  - oj(uj  - uj.,))  ^ (ioa) 


- 

u* 

U^  U‘)  + 

(n)  «x 

(“j)  ^ 

+ 

(fc) 

Dj  (uj,.,  • 

. 2 uj  * 

Uj-l) 

(10b) 

dt 

■ 

Uj  , . 

■ ''x  (“ 

' “j) 

* (fc)’  ^ 

, (“j)  ^ (“j) 

(10c) 

dUj 

3t  “ 

(K/a)  - 2 

Uj  * Uj. 

,)/(4x)X  ■ 

■ © 

(lOd) 

dU? 

iL  s 

dt 

- P2 

Uj  Uj  . P, 

("  - “j 

“j) 

(lOe) 

for  j * 1,  2,  ...  M where  M is  the  number  of  mesh  points  used  to  partition  the 
x-interval  of  interest.  The  existence  of  a solution  has  been  given  in  various 
books  for  systems  of  ordinary  differential  equations,  namely  Hille  (ref.  13), 
Caratheordory  (ref.  14).  The  most  general  conditions  for  existence  are  given  in 
reference  14;  but,  under  his  assumptions,  uniqueness  cannot  be  assured.  Stronger 

14 
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assumptions  are  necessary  and  these  are  Investigated  in  reference  13.  The 
equations  below  provide  a very  compact  method  for  rewriting  equations  (10). 


* f i(x,  t,  U) 

(Ha) 

^ u) 

(lib) 

+ f 3(X,  t,  U) 

(11c) 

= f 4(x,  t,  U) 

(lid) 

- f 5(x,  t,  U) 

(lie) 

with 

U'^(x,  0)  = Uj.  (x) 
which  may  be  rewritten  as 

_ U'  = f(x,  t,  U)  (12a) 

U(0)  * U"  (12b) 

Caratheordory  states  the  following  theorem  (ref.  14): 

If  the  n functions  f^  are  continuous  in  a closed  and  bounded  region  5"  of  an 
n + 1 dimensional  space  through  each  interior  point  (t°,  u°)  of  G,  there  exists 
at  least  one  continuously  differentiable  curve  * X^(t)  which  is  defined  on 
|t  - t°|  <a  and  satisfies  equation  (12).  Since  this  theory  says  nothing  about 
uniqueness  nor  how  the  solutions  depend  on  the  initial  conditions  nor  how  the 
solutions  vary  when  the  functions  f^  are  perturbed,  then  one  usually  requires 
that  the  function  f^  have  more  stringent  conditions.  These  fall  into  two  classes. 
The  first  is  concerned  with  bounding  the  partial  derivatives  of  the  f^  with 
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respect  to  the  and  the  second  is  Lipschitz  continuity  of  the  function  f^  in 
the  U variables.  Since  f is  a vector  function  in  one  defines  its  deriva- 

tive as  the  linear  transformation,  say  B that  satisfies  (ref.  15) 

Ltm  ft”)  - • 0 

IlhlH  IIMI 

, if  the  limit  exists  and  such  a transformation  B can  be  found.  The  symbolization 
h 0 implies  for  c > 0,  | |h| | <0  where  h c R".  The  matrix  B which  is  the 
Frechet  derivative  of  f is  related  to  its  Jacobian  . A theorem  due  to  Kamke 

(1930)  and  to  the  authors  Coddington  and  Levinson  (1955)  is  as  follows  (see 
reference  13): 

Let  f(x,y)  be  continuous  in  the  |x  - x# | < a,  | |y  - y| | jc  b and  satisfy 
||F(x,  z)  - F(x,  z)||  < w (lx  - Xgl,  ||y  -yaH)  where  W satisfies  the  following 


1. 

Let  H(r)  * j 

|h  eC'^  [0, 

t]  |h(0)  * 0,  and  Lim  = o[ 

ti-O  ^ 1 

2. 

W(x,  y)  > 0 

3. 

If  yi  > y2 

then  for  x 

g fixed,  W(Xg,  yj  > W(xg,  yj 

4. 

3c  > 0,  such  that  W(x 

, cx)  i c,  x£(0,  aj 

5. 

For  h EH(r) 

the  Lim  W(t,  h (t))  = 0 

then  the  only  solution  in  H(a)  of  z'  * W(t,  z(t))  is  z(t)  = 0 and  y'  « F(x,  y) 
has  a unique  solution  passing  through  (xg,  y,). 

Another  theorem  is  given  in  Hille  (ref.  13)  whereby  the  Lipschitz  condition 
on  F is  replaced  by  a condition  on  its  Frechet  derivation,  that  its  Frechet 
derivative  exist,  be  continuous,  and  be  bounded  in  the  domain  0.  The  full  theorem 
(ref.  13)  states  that  in  the  differential  equation  y'  ■ f(x,  y),  f:  R‘  x r”  r”, 
where  f is  continuous,  bounded  in  D * |(^,  y)_e_R^  1^  " *» 

My  ■ yoM  l 1 1^(^»  y)M  ^ M,  for  a11  (x,  y)  eO,  the  Frechet  derivative 

exists,  is  bounded  by  B',  for  small  values  of  | |y  - y,||  we  have  f(x,  y)  - 

(f,  yg)  ■ (y  - y#)  + o(||y  - yg||)  then  3!  solution  of  y'  « f(x,  y)  passing 

through  (xg,  yg)c  D. 
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One  may  obtain  a solution  to  system  (10)  by  assuming  those  parts  of  the 
system  due  to  "diffusion"  to  be  zero.  Under  these  circumstances,  the  results 
are  as  follows: 


* 9(x»  t,  U) 


Mi 

•C2 


Mi 

^3 


- Pj  U*  U*  + PjU®  (N  - U®) 

with  U**  being  constant  throughout  the  medium.  The  parameters  Pj  are  then 
determined  from  the  energies  of  the  electrons  and  holes  where  they  are  assumed 
to  be  a Fermi-Dirac  gas.  Under  this  formulation,  Kittel  gives  the  probabilities 
for  trapping  as  follows  (ref.  11): 

Pj.  » 1/^1  + exp  ({E  - £j')/(2  kbU-))) 

where 


dMi. 

at^ 


Ej'  * constant  (Uj"*)^^^ 

with  Sj  * constant  (Uj)^'^^  These  equations  were  modeled  under  a system  of  pro- 
grams called  Extended  Sceotre  (ref.  16).  These  routines  allow  one  to  solve 
systems  of  stiff  ordinary  differential  equations  as  well  as  electronic  circuits. 
With  only  a few  cards,  one  can  implement  the  above  system  without  having  to  write 
the  integration  routines.  The  implementation  used  is  given  in  the  appendix. 
Control  of  the  time  increment  is  by  SCEPTRE  and  before  going  to  the  next  time 
step  the  program  Iterates  the  solution  to  obtain  a "better"  solution.  If  after 
a fixed  number  of  tries  the  program  finds  nonconvergence,  then  the  program 
reduces  the  time  step  size  and  redoes  the  problem  with  best  known  data.  This 
type  of  control  is  necessary  when  stiff  differential  systems  are  encountered  to 
reduce  truncation  error  as  well  as  round  off  error  (by  reducing  the  number  of 
step  calculations).  The  system  also  calculates  the  Jacobian  matrix  noted  above 
and  notes  whether  the  matrix  is  of  full  rank  or  singular. 
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SECTION  V 
CONCLUSIONS 

There  are  two  problems  In  solving  systems  of  partial  differential  equations 
with  finite  differences:  consistency  and  stability.  The  system  of  equations 
(4)  or  (6)  exhibited  a problem  called  instability  which  can  radically  affect  the 
convergence  to  a solution.  In  using  a purely  explicit  formulation,  one  has  to 
use  such  small  At's  that  the  resulting  answers  for  times  corresponding  to  1000 
initial  At's  are  almost  entirely  error  because  of  round  off  error.  This  problem 
was  also  apparent  in  the  purely  implicit  formulation  because  the  answers  would 
still  not  converge.  The  answers  at  the  second  and  third  steps  were  grossly  in 
error  in  both  cases  (the  electric  field  was  given  as  10*®  volts  per  meter). 

Thus  there  was  an  error  in  the  modeling  of  the  problem  and  the  equations  were 
not  consistent,  hence  convergence  would  be  impossible. 

The  second  technique  exhibited  a different  problem  in  that  SCEPTRE  said  the 
matrix  B was  singular.  While  it  calculated  answers  for  the  first  time  step,  no 
answers  were  calculated  for  the  second  time  step.  Very  small  time  steps  were 
used  so  that  this  would  not  happen,  but  no  results  were  forthcoming. 

It  appears  that  the  linearization  used  in  the  second  technique  does  not  yield 
a useful  approximation  to  the  original  problem.  In  fact,  no  approximation  is 
found  at  all.  A more  exact  approximation  to  the  original  equations  would  perhaps 
change  the  singularity  of  the  Jacobian  matrix  and  thus  yield  a useful  approxima- 
ti on_  to  the _sqi utipri" of Jsystem  ( 6 ) . _ 
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APPENDIX 


LINEAR  APPOXIMATION  BY  ORDINARY  DIFFERENTIAL  EQUATIONS 


CIRCUIT  OCSC"^IRTn'r 

PYla  ELPCTUTt  PltUI''*  PY>=  fLtCrxON  CONCi-' NT‘*  a f luN . PYJ=  HOLt 
CONCCNTPar  (HM « PY<»=  Tt  ■'PERA  TUW^  t Rv;3  .•■iiimmc.R  OF  Trap  n PYbs 
internal  EMKKUY.  pr,=  HIi-hT  HANI*  ilOP  !►*.  EDRClNit  (■■iJNCriON. 

THIS  Formulation  assumes  constant  tcmpepatuk  and  ferhi  oas-type 

LAh  F.jP  CALOliLATINf-  CAPTUHE  PROHAm  IL  I T lES  . 

defined  Parametfos 

INITIALIZE  the  constants 
PL=  lO.E-o 
PEPS=  3.3a03?E-U 
PO=  l.bU?E-lR 
PNT»  I.IE23 
PK=  1 .3*li)62E-2Y 
PCs  5.AAlS<)E-3a 
PEN*0 
PT?sl .e-H 

PT3  = l.f-<5 

INITIALIZE  THE  VARIAbLES 
PYls  I.E6 

PYa=l.E23 

PY3»1.E2? 

PY5*  .izbSTnei? 

PYSsO 

SET  UP  THE  COEFFICIENTS 
PEFs  XEF (PC* (PY2»PY2) •• ( i . /D. > I 
PEFPs  xEFP  (PC* (PY3*PY3) ( I ./3. I ) 

PYAs  01 (PEPS.OYl tDPYl ) 

PM2S  am2  ( (PYA»»(-2.5)  l•(l.A36EA)*EXP(-^.SS5E‘•/AdS(PYl)  ) » 

PM3s  XM3(  (PYA**(-P.5)  ) •(B.SEA) *EXP(-A.7AEA/ASS(PY  n ) t 

PNIs  ANi ( (3. 1 S23t23) • (PYR**1 .S)*  FXP ( -PO* ( PEF ♦PEFP ) / ( 2 . ‘PK ‘PYm » I ) 

PGs  XG (PQ* (PnI»PNI-PY2»PY3) / ( ( PY2 .PY3 ♦PN I ♦PN I ) •PT2) ) 

PFE*  XFE ( l./t I .♦EXP  1 (PEN-PEF) / (PKOPYA)  ) ) ) 

PFPs  XFP(1./(I.*EXP( (pEn-PEFP) / (Pk*PY4  » ) ) ) 
define  the  nlFFEHENTIAL  EOUATIONS 
OPyI*  Xl (-“0*(P«2*PY2.pm3*Py3) »Py1/P£PS) 

0PY2X  x2(PG-(Py2/PT2) ) 

DPy3s  x3(PG-PY3/PT3) 

UPYS*  x5(-PFF»PYS*UPY2»PFP*0PY3*(PnT-PY5) ) 

0PY6*  xe>(PY6-PL*PYl*DpYl  I 
FUNCTIONS 

QKA.b.Os  (FDM(A,atC)> 

OUTPUTS 

PT I .Pr2.PY3.PYA .PYS.PYS.PEF.PEFP.PNI .PI  OT 
XSTPSZ.UPYl  .nPY2.UPY3.0PTS.0PY/S 
RUNCONTHOlS 

maximum  PPIrjT  points*  2000 
computer  time  lIHIT  * 20 
starting  step  size*  I.E-RO 

MINIMUM  STEP  SIZE*  l.E-RS 
STOP  time*  «0.E-9 
INTEGRATION  OOIITINE  * <PO 
RRITE  sImijlS  OaTA 
Subprogram 

function  fun  ( r.B.C) 

COMMON  /CNTRi.S/TImE 
X*  A*d*C 

if  (TI.mE.EO.O'  X*293.1«I 

FUn»X 

keTuRn 

ENO 
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B = B 
D ='D  = 

E = U®  = T 
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E' 

H = It 


Jn.  Jn 
n ’ p 


K 

K. 


M 
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^ ''j  “j  ^ 


R 

a 

f 


n(U^),  p(U®) 

q 


|A(x.  t)| 
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^^o 

U2»  ^3 
P 

'!'2.  Tj 

M II 
I I 

O 

3 
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ABBREVIATIONS  AND  SYMBOLS 


magnetic  induction 

displacement  field 

electric  fields 

total  energy  of  the  system 

energy  state  of  an  electron  or  hole 

magnetic  field 

electron,  hole  current  densities 
thermal  conductivity 
Boltzman's  constant 


number  of  mesh  points  partitioning  an  interval  [a,  b]CR 
total  number  of  traps 

probability  that  an  electron,  hole,  will  be  captured  by  a trap 

the  real  line 
thermal  diffusivity 
frequency 

electron,  hole  densities 
electric  charge 

last  upper  bound  of  t for  all  x the  closure  of  the  set  G 
time 

curve  in  space 

first  order  difference  in  x of  the  function  f = f(x)  - f(x  - Ax) 

is  an  element  of  permittivity 
wavelength 

permeability  of  free  space 
electron,  hole  mobilities 
total  electric  charge 
electron,  hole  average  lifetimes 
norm  (depends  on  the  space) 

Euclidean  norm 

impl ies 

there  exists 

there  exists  a unique 

containment  (eg.  AC  b for  x eA->x  cB) 
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